Dependence of exciton transition energy of single-walled carbon nanotubes on 

surrounding dielectric materials 
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We theoretically investigate the dependence of exciton transition energies on dielectric constant 
of surrounding materials. We make a simple model for the relation between dielectric constant of 
environment and a static dielectric constant describing the effects of electrons in core states, a bonds 
and surrounding materials. Although the model is very simple, calculated results well reproduce 
experimental transition energy dependence on dielectric constant of various surrounding materials. 
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I. INTRODUCTION 

Photoluminescence (PL) of single-walled carbon nan- 
otubes (SWNTs) has been intensively studied for elu- 
cidating their unusual optical and electronic properties 
due to one dimensionality^^^^^^^^. 
Since both of electron-electron repulsion and electron- 
hole binding energies for SWNTs are considerably large 
compared with those for conventional three-dimensional 
materials, the Coulomb interactions between electron- 
electron and electron-hole play an important role in 
optical transition of SWNTs 1 - 7 .^^^^^. Optical 
transition energies of SWNTs are strongly affected 
by the change of environment around SWNTs such 
as bundling 2 ^, surfactant suspension 7 - : 14 ' 24 and DNA 
wrapping 25 . Lefebvre et al. 7 reported that the tran- 
sition energies for suspended SWNTs between two pil- 
lars fabricated by the MEMS technique are blue-shifted 
relative to the transition energies for micelle-suspended 
SWNTs. Ohno et al^ have compared the PL of sus- 
pended SWNTs directly grown on a grated quartz sub- 
strate using alcohol CVD technique^ with SDS-wrapped 
SWNTs 2 . The energy differences between air-suspended 
and SDS-wrapped SWNTs depend on (n, m) and type of 
SWNTs [type I ((2n+m) mod 3 = 1) or type II ((2n+m) 
mod 3 = 2)2&2L2S]. 

Recently, Ohno et al. studied En transition energies of 
SWNTs in various surrounding materials with different 
dielectric constant, K cnv 29 . Observed dependence of En 
on K cnv for a (n, m) nanotube showed a tendency that 
can be roughly expressed as 



En — + ^m K env 



(1) 



where En denotes a transition energy when 

^env is in- 
finity, is the maximum value of an energy change 



of En by K onv , and a is a fitting coefficient in the or- 
der of 1, respectively. At this stage, the reason why the 
experimental curve follows Eq.(TT]) is not clear. 

In the previous theoretical studies of excitonic transi- 
tion energies for SWNT a 17 ' 19 : 21 ' 22 , a screening effect of a 
surrounding material is mainly described using a static 
dielectric constant k. However, since k consists of both 
K cnv and screening effect by nanotube itself, Ktube, ex- 
perimental dependence of transition energies on dielectric 
constants of environment can not directly compared with 
calculation a 17 : 19 ' 21 : 22 using the static dielectric constant 
k. In this study, we make a simple model for the relation 
between K onv and k. The calculated results of excitons 
for different K env reproduced well the experimental tran- 
sition energy dependence on dielectric constant of various 
surrounding materials. 



II. THEORETICAL METHOD 

A. Exciton transition energy 

Within the extended tight-binding mode l 21 : 22 ' 28 , we 
calculated transition energies from the ground state to 
the first bright exciton state by solving the Bethe- 
Salpeter equation, 



[£(k c )-£(k v )]<5(k c ,k c )<5(k v ,k v ) 

fiv(kX,k c k v )} *"(k c k v ) = a n ^ n (kX) 



(2) 
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where k c and k v denote wave vectors of the conduc- 
tion and valence energy bands and E(k c ) and E(k v ) 
are the quasi-electron and quasi-hole energies, respec- 
tively. f2„ is the energy of the n-th excitation of the 
exciton (n = 0, 1, 2, • • • ), and ^"(kckv) are the excitonic 
wavef unctions. The kernel K(k c k v , k c k v ) describes the 
Coulomb interaction between an electron and a hole. De- 
tails of the exciton calculation procedure is the same as 
presented in Ref a 21 ' 22 ' 28 . 
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The exciton wavefunction |\&™ > with a center-of-mass 



momentum q(= k c — k„) can be expressed as 



I*™ >= 



k 



■^kc,(k-q)u C kc C (k-q)w|0 >, 



(3) 



where ZV 



J kc (k-q)« IS the eigenvector of the n-th (n = 
0, 1, 2, • • • ) state of the Bethe-Salpeter equation, and |0) 
is the ground state. Due to momentum conservation, the 
photon-excited exciton is an exciton with q « for par- 
allel excitations to the nanotube axis. In this Letter, we 
calculate the n = state of q = exciton for each (n, m) 
SWNT. 



B. Dielectric screening effect 

In our calculation, the unscreened Coulomb potential 
V between carbon n orbitals is modeled by the Ohno 
potentiate. We consider the dielectric screening effect 
within the random phase approximation (RPA). In the 
RPA, the static screened Coulomb interaction W is ex- 
pressed a&i 7 - 



W = Vy«e(q), 



(4) 



where e(q) is the dielectric function describing effects of 
the polarization of the n bands. K is a static dielec- 
tric constant describing the effects of electrons in core 
states, a bonds, and surrounding materials. In the cal- 
culation, we directly calculate only the polarization for 
the 7r band, and the effects of electrons in core states, 
a bands, and surrounding materials are represented by 
a single constant k. In the most accurate expression, 
the inhomogeneous and nonlocal dielectric response of 
the nanotube itself and the surrounding materials should 
be considered. However, it is not easy within extended 
tight binding method. In this study, instead of treating 
the complicated dielectric response including surround- 
ing materials, we make a simple model for a relation be- 
tween the static dielectric constant k and K cnv to obtain 
the En dependence on n cnv . 



C. Relationship between k and K en v 

Figure[T]shows a schematic view for the model relation- 
ship between n and K env . Here we consider the screening 
effect related to linear combination of the screening 

of nanotube itself and the surrounding material 



1 

K 



Ctu 



be 



Co 



K tubc 



(5) 



where Ktubc is the dielectric constant within a nanotube 
except for the 7r bands, and Ctube and C cnv are coef- 
ficients for the inside and outside of a nanotube, re- 
spectively. As shown in Eq.([T]), the transition energies 
observed in the experiment 29 indicate that there is a 





FIG. 1: Schematic of the connection of the net dielectric con- 
stant k and the dielectric constant of the surrounding material 
«env and the nanotube itself fttube 



limit valued when 



oo. Hence, when K e 



°o, C cnv /n cnv can be removed from Eq.©, and 1/k is 
expressed by the limit value as 



1 

K 



Ctube 



/^tubc 



1 



i (t e 



oo) 



(6) 



^tubc 



where 



is the limit value of the net dielectric con- 
stant k when Konv is infinity. Since electric flux lines 
through inside of the nanotube remain even when K env 
— > oo, we assume there is a certain value of k (ftj^ bo ) 
that corresponds to the situation when dielectric screen- 
ing by surrounding material is perfect and only dielectric 
response of the nanotube itself contributes to the net 
screening effect. 

Replacing C tu bc/Ktube by K^ bc , Eq.© is modified as 



1 



1 



K ^tubc 



(7) 



Next, we imagine that the SWNT is placed in the vac- 
uum, which corresponds to n = K vac and K on v = 1, and 
then Cenv can be expressed as 



C'e 



1 

-vac 



l 



"tube 



(8) 



where K vac is the static dielectric constant not for the vac- 
uum, but for the situation that the nanotube is placed 
in the vacuum. We now express k as a function of K env 
through two parameters K^ he and K vac , whose values can 
be estimated from the following discussions. In the pre- 
vious paper o 17 ! 19 ! 21 ! 22 , k value is put around 2 to obtain a 
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FIG. 2: (a) The En energy for a (9, 8) SWNT as a function 
of k. (b) SEn dependence on K cnv . Inset in (a) shows the En 
dependence up to k = 100. In (b), circles denote the experi- 
mental data and solid curves denote the calculated results of 
Eq.Hnj for K vac = 1.0 (black), 1.5 (red) and 2.0 (blue). 



good fit with experiments for SWNTs with surrounding 
materials. Jiang et al£^ have compared the calculated 
results with the results for the two photon absorption 
experiments^, and obtained the best fit using k — 2.22 
for SWNTs in a polymer matrix. Here, since K vac is for 
nanotubes without surrounding materials, K vac should be 
less than about 2 and close to 1 due to vacancy of inside 
of the tubes. With regard to K^ be , according to the ex- 
perimental results^ ' i 29 , transition energy change due to 
change of surrounding materials is at most 30-100meV. 
FiglUa) shows the calculated En energy dependence on 
k for a (9,8) SWNT in a small k region, while the inset 
shows the En dependence up to K onv = 100. As shown 
in FigJSJa), variation of k that yields the transition en- 
ergy change of 30 to 100 meV is about 1 to 3 when k is 
around 2. Therefore, the value of K^J bo should be around 
2 to 3 and that of K vac should be around 1 to 2. 



D. Dependence of excitation energy on K en v 

As shown in Fig^a), the calculated E\i energies de- 
crease with increasing k. This is mainly due to the fact 
that the self energy (e-e repulsion) always exceeds to 
the e-h binding energy and that the both interactions 
(e-e and e-h) decrease with increasing n. The En al- 
most linearly depends on k around the small k region. 
We checked that the linear dependence is universal for 
all (n, m)'s for diameters more than 0.7 nm. Assuming 
the linear dependence, variation of the excitation energy 
SEn = En — En(K env = 1) for the small k region is 
approximated by 



where A nm is the gradient of SEn near the small k region 
for each (n, m) type. After we transform k using the 
relationship of Eq. ([5]) , Eq. ([9]) is modified as 



SEn = -A nm {nZ\ 



be 



=)(- 



1 



+ («t~ 



be 



(10) 

A nm (n^ hc — K vac ) corresponds to the maximum value 
of SEn when K onv — > oo, which corresponds to the value 
of coefficient in the fitting curve of Eq.{T]). For 

(9, 8) SWNT, the fitted value to the calculated results 
for A nm is 33 meV and A^ obtained by the fit to the 
experiment^ 9 - using Eq.([T]) is 36 meV, and K^ hc — K vac 
should be around 1. The values for n^ hc and K vac 
are consistent with the values conventionally used for 
SWNTs in dielectric materials!^! 2 ^. 



III. RESULTS AND DISCUSSION 

Figure^b) compares SEn for a (9,8) SWNT depend- 
ing on K cnv by the experiment (solid circles) and the cal- 
culated results (lines) for n vac = 1, 1.5, 2.0 using Eq.^U]). 
As shown in FigJ^b), the qualitative shape of theoretical 
curves are in good agreement with the experiment and 
not affected so much by the change of K vac . Since the ex- 
act value of K vac is unknown, we hereafter set n vac = 1.5 
for each (n, m) SWNT. For the (9, 8) SWNT in FigOJb), 
K tubc = 2-7 and K vac = 1.5 are fitting values. These 
values are consistent with the discussion in the previous 
section. After setting K vac = 1.5, Eq. lfTU)) turns to be 



SEn 



A nm (^e 



i) 



>W« bo " * vac ) + 1/1-5 



(11) 



Thus, we express SEn as a function of « env with one 
parameter (K^ be — K vac ). 

Figure[3]^a) shows the calculated values of A nm for each 
(n, m)'s. Family pattern of (2n + m = const.) family 
is drawn with the 2n + m values by dotted lines. We 
found a slight diameter dependence and relatively large 
chiral angle dependence of A nm for type II SWNTs (blue) 
compared with type I SWNTs (red). The type II SWNTs 
with larger chiral angles tend to have larger value of A nm . 
For a convenient use of Eq. (flQ|) . we give a fitting function 
of A nm meV as 



A nm = A + Bdt + (C + D/d t ) cos 3(9, 



(12) 



SEn 



A nm (k, k ), 



(9) 



which gives the average (maximum) error of 
±2meV(8meV) for type I, and ±2meV(5meV) for 
type II SWNTs. The fit curve is shown in Figi^a) by 
solid lines. Here dt (nm) is the diameter of nanotube 
and 9 is the chiral angle 2 ^. The values of (A, B, C, D) 
are (36, -4, 0, 0) and (33, -3, 6, 7) for type I and for type 
II SWNTs, respectively. 

In order to expand our result to many (n, m) SWNTs, 
we need a function to describe («^ be — K vac ). It is impor- 
tant to note that K^ bo should depend on the diameter. 
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FIG. 3: (a) Calculated values of A nm for each (n, m) SWNT. 
Open (red) and solid (blue) circles correspond to type I and 
type II SWNTs, respectively. Solid lines denote the fit curve 
by Eq. JT2J . (b) «J£ be - K vac vs l/d 2 t . The values of n^ be - K vac 
are obtained by the fit of Eg. (|10fl to the experimental data 
for each (n, m). 



An exact function should be calculated by taking into ac- 
count the Coulomb interaction considering induced sur- 
face charge at the boundary of the nanotube and sur- 
rounding materal for an e-e or e-h pair for each (n, m) 
SWNT. Instead of calculating this complicated function, 
here we roughly estimate the (K^J bc — K vac ) as a simple 
function of diameter d t , since (Kj^ bc — K vac ) should depend 
on the cross section of a SWNT. As shown in Fig[3Jb), 
( K tube — KVa °) is roughly proportional to , 



''tube 



E 



with the coefficient E = 1.5±0.3 nm 2 . Here («:^ bc 



(13) 

K VaC ) 

is obtained by the fit using Eq. (fT0)) and A nm calculated 
for each chirality. FiglH^b) clearly shows that our cal- 
culated A nm well describes the chiral angle dependence 
of 8 En and that the remaining diameter dependence is 
understood by («££ be - n vac ) through 1/d 2 . This 1/df 
dependence implies that K^ be depends on the volume of 
inner space of the nanotube. Although the number of 
experimental data available for the fit is small and se- 
lection of this function is arbitrary to some extent, it is 
reasonable that 1 / K^ bc increase with the increase of the 
diameter, since l/K^ bo corresponds to the Coulomb in- 
teraction through the inner space of the nanotube. In 
order to find an accurate form of the function, future ex- 
periments and theoretical studies are definitely needed. 

Figure[l]shows SEn as a function of K cnv for (a) the ex- 
periment and (b) the calculation using Eq. (fTTj) and (fl"3)) . 
FiglUc) compares 8Eu for the experiment and that for 
the calculation with the same K env values. The same sym- 
bols for an (n, m) are used in three figures of FigSI De- 
tails of experimental data will be published elsewhere^. 
Although our treatment is very simple, the calculated 




FIG. 4: The transition energy dependence plotted as a func- 
tion of K onv . (a) experiment and (b) calculated results are 
indicated by (a) symbols and (b) solid curves. In (b), (n,m) 
for each curve is indicated by a symbol on the curve, (c) 
Comparison of SEn for the experiment (SEn (experiment)) 
and calculation (SEn (theory)). A dotted line indicates the 
line of SEn (experiment) = 5En (theory). Open (red) and 
solid (blue) symbols correspond to type I and type II SWNTs, 
respectively. The data in the dotted circle are the data for 
k ci1v = 1.9 2 — (see text). 



curves for various (n, m) SWNTs well reproduce the ex- 
perimentally observed tendency for each (n, m) SWNT, 
and the degree of difference between each (n, m) type is 
also in good agreement with the experiment. As shown 
in FigQJc), 5E\\ (theory) is in a good agreement with 
5En (experiment) except for several points indicated by 
a dotted circle in the figure, which correspond to a case 
for the smallest /c cnv = 1.9 (hexane) except for K onv = 1 
(air) in the experimental data™. The value of K cnv = 1.9 
for hexane is adopted as the dielectric constant for the 
material, in which the dipole moments of liquid hexane 
are not aligned perfectly even in the presence of the elec- 
tric field. Since K cnv = 1.9 is a macroscopic value, a local 
dielectric response might be different from the averaged 
macroscopic response. If the local dielectric constant 
near SWNTs becomes large (for example, K onv w 3), the 
fitting of FigQJc) becomes better. We expect that the 
dipole moments of a dielectric material might be aligned 
locally for a strong electric field near an exciton, which 
makes the local dielectric constant relatively large. This 
will be an interesting subject for exciton PL physics. 
Since the difference of A nm between each (n, m) type de- 
creases with increasing the diameter, it is predicted that 
the amount of variation due to the change of K C nv mostly 
depend on diameter in the larger diameter range. Thus a 
PL experiment for nanotubes with large diameters would 
be desirable for a further comparison. 



IV. SUMMARY 

In summary, the dependence of exciton transition en- 
ergies on dielectric constant of surrounding materials are 
investigated. We proposed a model for the relation be- 
tween dielectric constant of the environment and a static 
dielectric constant k in the calculation. Although the 
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model is quite simple, calculated results well reproduce 
the feature of experimentally observed transition energy 
dependence on dielectric constant of various surrounding 
materials, and various dt and 9. 
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